UCSD/PTH 02-17 
INT-PUB 02-43 
hep-ph/0208034 

Testing Factorization in B D^'>X Decays 

C.W. Bauer/ B. Grinstein,^ D. Pirjol/ and I.W. Stewart^ 

^Department of Physics, University of California at San Diego, La Jolla, CA 92093* 
^Institute for Nuclear Theory, University of Washington, Seattle, WA 98195^ 

Abstract 

In QCD the amphtude for D^*^^tt^ factorizes in the large Nc limit or in the large energy 

limit Q ^ Aqcd where Q = {mh,nic,mb — me}. Data also suggests factorization in the processes 
B D*7r^7r^7r^7r'^ and B D*u}it~ , however by themselves neither large N^ nor large Q can 
account for this. Noting that the condition for large energy release in B^ D^tt^ is enforced by 
the SV limit, nih S> mf, — rric ^ A, we propose that the combined large N^ and SV limits justify 
factorization in i3 — > D^*'^X. This combined limit is tested with the B — > D*X inclusive decay 
spectrum measured by CLEO. We also give exact large Nc relations among isospin amplitudes 
for B and B D^*W^*^X, which can be used to test factorization through exclusive 

or inclusive measurements. Predictions for the modes B D^*'7nr, B D(*^KK and B 
are discussed using available data. 
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I. INTRODUCTION 



There are two limits of QCD in which factorization in —>■ D^tt^ (and related) decays 
can be proven rigorously. In the large limit, in which one takes the limit of infinite 
number of colors, one can show that all nonfactorizable contributions are suppressed by 
In the large energy limit, Q ^ Aqcd? where Q denotes any one of the large scales 
rrib, rric or E^^, factorization has been proven rigorously []1| using a soft-coUinear effective 
theory . In this case the corrections to factorization are suppressed by powers of 1/Q . 

The large Q proof of factorization has its roots in a non-rigorous argument by Bjorken 
that justifies factorization using color transparency. Here the requirement is that the pion 
velocity in the B rest frame be ultrarelativistic, so that the small color dipole of quarks that 
forms the vr" grows to hadronic size A only after a time delayed hadronization time 7A much 
larger than the hadronic size of the B meson. Although not systematic, color transparency 
therefore seems to suggest that violations to factorization are order I/7. While both large 
Nc and large Q lead to factorization in a particular limit of QCD, it is not clear which 
is the most important for the phenomenologically observed factorization in real QCD, in 
which neither large Nc, large Q, nor large 7 are exactly realized. Each method gives similar 
predictions for the decay B^ —>■ D^n^ and, while none predict the decay B^ D^ir^, they 
all render it suppressed. Hence these decays are not very useful for distinguishing among 
these arguments for factorization. 

If large Nc were the only requirement for factorization, one would expect factorization 
to hold equally well in D decays and K decays. Factorization does not work very well in 
D decays 0, Q and is even more strongly violated in K decays. Clearly large Nc by itself 
does not explain phenomenological factorization. A great deal of effort has been devoted to 
explaining this aparent puzzle |^ ||. 

On the other hand, if the large energy limit was solely responsible for explaining fac- 
torization in B decays, one would expect corrections to grow as mx/Q, where mx is 
the hadronic mass produced by the {du)v-A current. Recently, an analysis of the decays 
B D* TT'^ 7T~ 71^ and B D*ujtt~ was performed, investigating the applicability of fac- 
torization as a function of the invariant mass of the light hadrons 0. Using information 
from r decays it was shown that up to mx ~ 1.7 GeV there is no indication of violations 
to factorization, which indicates that the large Q limit can not be solely responsible for 
explaining factorization. However, this data can not be explained solely from the large Nc 
limit of QCD either. In order to calculate the predictions from factorization, Ref. used 
the factorized form 

{XD^*^\{ch)v-A{du)y^A\B) = {D^*^\{ch)v-A\B){X\{du)v-A\Q) . (1) 

This neglects the possibility that part of the hadronic state containing the light particles 
can be created from the b ^ c current. In fact, large Nc only gives 

{XD^*^\{cb)v-A{du)v-A\B) = Yl {D^*^X'\{cb)v-A\B){X"\{du)v-A\0) , (2) 

X',X" 

with X' and X" adding to give the final light hadron state X. The authors of Ref. 
addressed this issue by arguing that the contributions from non-zero X' may be numerically 
small and that for B —>■ D*ujn~ this can be tested using data from semileptonic B — > Duiu 
decays. These effects can also be tested in the decay B^ —>■ D*^TT~^7T^n^7T~ for which 
X' = {tt"*" , 7r^7r"*"7r~} and X" = {tt+tt^tt^, tt^}. The CLEO measurement [|iy] finds a small 
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but nonvanishing branching ratio for this mode Br{B'^ D*^7t^tt^tt tt ) = (0.30 ± 0.07 ± 
0.06)%. 

It is interesting to note that factorization in Eq. ([1|) is justified in the small velocity (SV) 
limit, rrih, rric ^ rrih — rric ^ Aqcd in|- shown in Ref. |T2[, in the SV limit the inclusive 



B — > Xciu branching ratio is saturated by D and D* final states. Therefore, in the SV limit 

{X'D^*^\icb)v-A\B) = {D^*^\icb)v-A\B)Sx',o- (3) 

Thus, Eq. reduces to Eq. (|1]) in the simultaneous limit of large Nc and SV and this 
combined limit is therefore capable of explaining the results of P]. 

Although more stringent than the large Q limit, the SV limit implies large Q and is 
therefore fully consistent with it and also gives additional predictive power. For B D^*"* 
transitions the relevant kinematic limits can be summarized as 



HQET 
Large Q 
SV limit 



nib, rric > Aqcd , 

rrib, rric, mb-rric':^ Aqcd , 

m^, rric nib - rric Aqcd , (4) 



each of which is a subset of the one above. The requirement that the light degrees of 
freedom in the D^*'> can be described by HQET requires that Ed ~ rric or equivalently that 
the B and D velocities have v ■ v' of order one. This implies that rrib — rric ~ \f^f^M^c or 
rrib — rric ^ y/^b^c which are allowed scalings in the large Q and SV limits respectively. 

The purpose of this paper is twofold. In section |I| we show that even the more generally 
factorized form of the amplitude in Eq. (H) leads to experimentally testable predictions 
which are distinct from those following from large Q factorization. The B — > D^*^Xu decays 
with charge eigenstates are parameterized by four independent isospin amplitudes. In the 
large Nc limit, they are given in terms of only two reduced matrix elements, compared with 
just one in large Q factorization. A similar result holds for B DDX decays, for which 
we prove a similar reduction in isospin amplitudes from 7 to 5, in contrast with two in the 
combined large Nc and SV factorization. In section |Tl| we consider the simultaneous limit of 
large Nc and SV and calculate the inclusive differential decay rate B ^ D*X in this limit. 



The results obtained are compared with data available from the CLEO collaboration [13 . 



II. LARGE Nc RELATIONS FOR B ^ DX 

The final state X in the decay B — DX may have charm number —1 or 0, depending 
on whether the underlying weak decay is 6 — > ccs or b ^ cud. We analyze these separately, 
correspondingly labeling the final states Xc and Xu- Taking into account the fact that the 
B and D states belong to isospin doublets and that the weak Hamiltonian 

Hw = -^VcbV*a\Ci{iJ,) {cu)v-Aidb)v-A + Ciifx) {cb)v-Aidu)v~A\ (5) 

transforms as an isotriplet with (/, /s) = (1,+1), one obtains that the final light state X^ 
can have isospin Ix^ = 0, 1,2. Specifying the isospin of the Xu system together with the 
total isospin of the final hadronic system gives rise to four isospin states: [DXo]i, [DXi]!, 
[DXi] ^ and [-DX2] 3 . The corresponding reduced isospin amplitudes are 

ai = {[DXo]i\\nw\\B) , 61 J = {[DXi]i^i\\nw\\B) , C3 = ([DX2] 3 ||Hvk||5) . (6) 
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The isospin amphtudes depend on X through its hadronic content and particle momenta. 
Squaring the amplitudes and summing over the isospin of the X„ state one finds for the 
rates with B and D charge eigenstates 

T{Bd ^ D^x-) = j: %,{x) + hy.ix) ' + j:Imx)\' 



X 5' 



^ 3 ' 3 

T{Bd - D'Xl) = E ?l«i/2(X)P + E H + h„{X) " + E ^|c3/2 WP 



3 ^9 ^ 

T{B- D'X-) = E |&3/2(X)f + E ^|C3/2(X)P 

X X ^ 

T{B- ^D^X--)=J:%V2{X)\\ (7) 

The phase space factors are implied. The corresponding sums over X in the rate formulas 
include summation over the final states and integration over the phase space. 

In the large Nc limit these amplitudes simplify considerably. We assign the usual Nc 
power counting to the Wilson coefficients Ci(mfe) ~ 1/iVc, C2{mf,) ~ 1, in agreement with 
their perturbative expansion in ag This gives that to leading order in 1/Nc the B DX 
amplitude factors as 

A^oiB ^ M„) = ^V,,V:,C2 E {DX'\cb\B){X"\du\0) + 0(l/iV,) (8) 

with Xu = X' + X". Since the state X" has to be in a state of isospin |/, I^) = |1, —1) and 
the state X' can only have isospin / = or 1, all possible decays B DX^ are determined 
in the large limit by two amplitudes Fi defined as 

Fo{X)= E ([M^]i||c&||5)(Xr||J«||0), (9) 

X'X" 



FAX) = E ([M;].||c6||5)(xri|rf«||o) 



X'X' 



The four isospin amplitudes (|^) can be written now in terms of Fq and Fi as^ 

ai/2 = -Fi, 6i/2 = Fo + y2Fi , 63/2 = ^0--^^!, C3/2 = y^^i • (10) 
Using these expressions in the rate formulas, one finds the large A^^^ predictions 



rNo(-B<i - 


- D^x;:) 


= T.\Fo{X) + ^F,{X)\' + - 


Fi{X)\' 


{Bd - 


. D^X^) 


= E2|Fi(X)p 

X 




rNo(5- - 


- D'X~ 


)=El^oW--^Fi(X)p + l 


Fi{X)\' 


rNo(5' - 


- D+X- 


-) = E2|Fi(x)r. 





(11) 



X 



^ Although in general Fo^i{X) ~ 0{N^), for special states X some of these amplitudes may vanish. For 
example, taking X = n gives i^i(7r) = 0. 
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The two independent amphtudes can be extracted from the combinations of rates 

r{Ba ^ D+X-) + r{B- ^ D^X-) = J2 2|i^i Wr + 2|Fo(X)|2 (12) 

X 

TiB, ^ = - D^X--) = 5:2|Fi(X)r . (13) 



X 



The corrections to the predictions in Eq. ( pil| ) come at order l/N^. and are parameterized 
by the subleading term in Eq. (^. This can be written again in a factorized form as 

A„^XB - M„) = ^K^KT. (Ci + ^) Y:jDX'\cu\Q){X"\dh\B) + ©(l/iV^) . (14) 

For this case both X' and X" can have isospin J = 0, 1. Thus four isospin amphtudes are 
required to describe all matrix elements (|Hp, which shows that in general no reduction in 
the number of isospin amplitudes persists beyond leading order in 1/Nc. 

The large energy limit can be applied to these decays if the final state X contains only 
particles which form a jet with energy much larger than its invariant mass (this includes 
the case of just a single light particle with large energy). In this case, the amplitude can be 
factored in a way similar to the large Nc limit as 

A^{B ^ DX^) = %V,,V:, (C2 + ^) {D\cb\B){X\du\0) + O L(g), ] (15) 



where mx is the invariant mass of the state X and n > 0. Since the isospin of X is 
constrained to be / = 1, the only nonvanishing isospin amplitudes in this limit are bi/2, 
This prediction can be tested for example by measuring the decay B~ D^tt~7i~ as a 
function of the angle between the two energetic pions, which should vanish as this angle 
decreases. 

The above analysis can be carried through in an identical manner for B D*X decays. 
In the combined large Nc and SV limit the amplitude {D^*^X\cb\B) vanishes for X 7^ 0. 
Thus, Fi{X) = for all final states X. This prediction is similar to the one from large 
energy factorization, but without the kinematic requirement on the light hadronic system 
X. Separate measurements of the 4 rates with B, D charge eigenstates would allow a test of 
these predictions, and distinguish between large factorization, large energy factorization 
or the combined limit of large and SV . In the next section we discuss a partial test along 
these lines, making use of the present limited experimental information on B ^ D^*^X 



available from CLEO 13 



We stress that the large Nc relations hold not only for the inclusive mode, but also for 
states X with fixed hadronic content. For example, taking X = tttt gives that in the large 
Nc limit, the amplitudes for the two pions in i? ^ D^*\TnT)i to be emitted in states of 
isospin 1 = and 2 are related. Eq. (p]3| ) gives a relation among rates to leading order in 
1/Nc (note that r{Bd D^*^'^n'^n^) ~ 0{1/N^)) 

T{Bd DW°7r+7r-) = T{B- D^*^^^-^-) + 0{1/Nc) . (16) 

The branching ratios of these modes have been recently reported by the BELLE Collab- 
oration For the D modes they are Br{Bd D^Ti+n-) = (7.5 ± 0.7 ± 1.5) x lO"'' 
and Br{B' D+n'Ti-) = (1.07 ± 0.04 ± 0.16) x 10-^ and for the D* modes Br{Bd 
D*\+n-) = (6.2 ± 1.2 ± 1.7) x 10"^ Br{B- D*+n-n-) = (1.24 ± 0.07 ± 0.22) x lO'^. 
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The data agrees fairly well with the large Nc prediction for the D modes and suggests 
larger 1/Nc corrections for the D* modes. This agreement permits a determination of 
|Fi(7r7r)p (integrated over phase space) using Eq. (p!3[). A subtraction can then be performed 
to extract |Fo(7r7r)p using (0). The widths on the left side of Eq. (|1^) are dominated by 
the two-body B Dp mode[|T^]. One finds (in units of Br) 

2|Fo(7r7r)p = (20.4 ±2.3) X 10"^ 2|Fi(7r7r)p = (9.1 ± 1.1) x 10"^ (17) 

The large observed enhancement of Fq over Fi is a consequence of the inequality among two- 
body modes T{Dp) ^ r(D**'7r). Factorization for these two body modes can be explained 
either in large N^. or in the large Q limit. The suppression of Fi over Fq is in agreement 
with the SV limit. 

Another test of the large Nc relations is obtained by taking X = KK, for which data is 



available from BELLE The isospin of the kaon pair can be only and 1, which gives 
C3/2{KK) = 0. The large Nc relations (|l3) imply Fi{KK) = 0, which shows that in this 
limit the B —>■ D{KK)j=o amplitude is suppressed, and the decay rates satisfy 

r(5d ^ D+K-R^) = T{B- D'^K-R'^) + 0{1/Nc) 

r{Bci D^K^R^) = r{Bd ^ D^K+K-) = 0{1/N^) . (18) 



The first prediction agrees well with the BELLE results []T6| Br{Bd D^K K*^) = 
(8.8 ± 1.1 ± 1.5) X 10~^ Br{B- D^R-R*^) = (7.5 ± L3 ± 1.1) x 10^^ and Br{Bd 
D*+K-R*^) = (12.9 ±2.2 ±2.5) x 10-^ Br{B- D*^K-R*°) = (15.3 ± 3.1 ± 2.9) x 10"^ 
Next we turn our attention to the b — > ccs process. For this case, large Q arguments can 
not be used to justify factorization, which leaves large Nc as the sole possible explanation. 
The Hamiltonian responsible for these decays is identical to (|^) with the substitution 
u ^ c. We neglect the penguin operators because of their small CKM factors and Wilson 
coefficients. We first consider the case in which experiments may tag the B meson and dis- 
tinguish the charm from anti-charm. There are two isospin amplitudes Hq, hi corresponding 
to the two possible values of the isospin I of the state Xc 

ho{X) = {[DXcoh_\\n'^\\B) , hi{X) = {[DXci]i\\n'^\\B) . (19) 
The corresponding rates are given by 

-j:\ho{X)\' + hh{X)\' 

X 

' = j:Imx)\' 

X 



r{Bd- 


- D^X;) 


r{B,- 


^ D'X!) = 


T{B- - 




T{B- - 


^D'X-) 



Y^\h,{X)\' + \\hi{X)\\ (20) 



X 



No simplifications are expected for these modes in the large Nc limit. Due to the identical 
isospin structure of their Hamiltonian, similar relations can be written down for the semi- 
inclusive semileptonic decays B DXeu, in terms of another two amplitudes gQ^i{X). 
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Table I. Isospin decomposition of the most general B DDXj decay amplitude into 
states Xj with isospin I. 



The definition of the SV hmit for decays containing two charm quarks in the final state 
is somewhat different from the one introduced in Eq. (|^). Requiring that both charmed 
hadrons move slowly gives 

Generalized SV limit: m^, rric^ mi, — 2mc ^ Aqcd ■ (21) 

In this combined large and SV limit, the state is produced by the (sc) current and 
must have isospin 0, which requires the amplitude hi{X) to vanish hi{X) 0. 

Making explicit the charm and anticharm in the final states gives many more modes. 
This type of analysis is necessary if the experimental inclusive measurement relies on the 
presence of a charmed meson and separates them only according to whether they are charged 
or neutral. Taking into account the fact that the Hamiltonian responsible for these decays 
is an isospin singlet / = 0, one finds seven independent isospin amplitudes describing these 
decays. Three reduced amplitudes describe decays into nonstrange D mesons 

ai = {D{DXi/2)o\\n'^\\B) a2 = {D{DXi/2)i\\n'^\\B) = {DiOXsMW^WB) (22) 

and another four reduced amplitudes parameterize decays into Dg 

h = {D{D;XoU\n'^\\B) 62 = {D{D;Xi)i\\n'^\\B) 

Cl = (D(D+Xo)o||H^||S) C2 = {D{DtX,),\\n'^\\B) . (23) 

Finally, another amplitude d describes B DfDj X1/2 decays. The isospin decomposition 
of the most general B DDX decay amplitude is shown in Table I. 

In the large Nc limit, the amplitudes for B — > DDX factor in a similar way as for 
B — > DXu (H), (|14D- Keeping terms up to 0{1/Nc), the factorizable terms read 

A{B^DDX) = ^VcbV:,C2 J2 {DX'\cb\B){DX"\sc\0) (24) 
v2 x',x" 
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Mode 


Decay rate 


Br X 10-3 


Bd 
B 


D+D^K- 




(1.7 ±0.3 ±0.3) 
(1.8 ±0.7 ±0.4) 


Bd 

B- 


-> D+D-R^ 




(0.8l[;|±0.3) 
(1.9 ±0.3 ±0.3) 


Bd 

B-^ 


^ D+D-R- 




(0.8 ±0.4 ±0.2) 
(0.0 ±0.3 ±0.1) 



Table II. Large Nc predictions for a few B DDK modes. The second column 
gives the large Nc prediction and the last column shows the branching fractions for 



the corresponding modes measured by the Babar Collaboration |17]. The consecutive 
pairs of modes have equal widths in the isospin limit. 



^K,V;: (Ci ± ^) j:jDDX'\cc\0){X"\sb\B) + ©(l/iV^) 



Taking into account isospin constraints on the intermediate states X',X", one finds that 
the 3 amplitudes ai_3 are given by 2 independent amplitudes Aq^i at 0{N^) (corresponding 
to X' in the first term having isospins / = 0, 1) and another 2 amplitudes Bq i at 0{1/Nc) 
(corresponding to X' in the second term having isospins / = 0, 1) 

ai = V2Ao - ^{Bo + B,) , = A^ - ^{B, - ^B^) , as = -V2A, + ^^^^ .(25) 

At leading order in 1/Nc there are 2 relations among amplitudes 

a,{X) = -V2a2{X) , c^{X) = -^C2(X) . (26) 

No such relations exist among the coefficients 6j. These predictions imply the large Nc rate 
relations (separately for each of the four channels B — > D^*^D^*^X and summed over the 
isospin of X) 

T{Bd D^D^X) = T{Bd D^D-X) = T{B~ D+D^X) = T{B~ X) 
T{Bd D+D^X) = T{Bd D+D-X) = T{B- D+D^X) = T{B- X) . 

(27) 

No data exists at present on any of these inclusive modes. It is possible to test the large 
Nc predictions on the exclusive modes B £)(*)£)(*)}{(*) using data recently available from 
Babar [1^. Taking X = K^*^ requires a^^K) = since the isospin 3/2 is not allowed. The 



large Nc relation (^) predicts a2{K) = 0{1/Nc), which implies the rate relations shown 
in Table II. Comparing these predictions from large Nc with the data one finds reasonable 
good agreement within the theoretical and experimental uncertainties. 

The large Nc predictions (p6D can be again contrasted with the more constraining results 
obtained in the combined large Nc and SV limits, according to which, among the isospin 
amplitudes ai,bi,Ci, d, the only nonvanishing ones are Oi and bi. For the modes in Table II 
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the SV hmit predicts a vanishing Bq. However, since this amphtude is subleading in 1/Nc, 
these modes can not be used to test the SV hmit. Therefore separate measurements of other 
modes in Table I are needed to distinguish which of these limits (if any) are actually realized 
in nature. 



III. FACTORIZATION FOR B D*X 

As a further test of factorization we propose studying the inclusive hadronic decay B — > 
D*X as a function of the invariant mass of the state X, = m^. CLEO has measured this 
decay spectrum from m"^ ~ out to the maximum hadronic mass = [rriB — mjy*)'^ c:^ 



10.7 GeV |T3[. Thus, the inclusive spectrum allows for a test of factorization over a much 
larger range than the ~ to 3.2 GeV^ region considered in Ref. 0. 

As explained in the introduction, the inclusive B —>■ D*X amplitude in the large Nc limit 
is given by 

MB D*X) = E CM{D*X'\{cb)v-A\B){X"\{du)v.A\0) , (28) 

x',x" 

with a sum over X = X' + X" in the rate, plus additional analogous contributions involving 
the {sc)v-A current and Cabibbo suppressed terms. Imposing the SV limit as well gives 
rise to two important simplifications to Eq. (PB]). First, as in Eq. (|^), the contributions 
from hadrons X' produced together with the D* from the {cb)v-A current are suppressed. 
Second, the contributions from non-prompt D*, that is, D* that arise from production of 
higher resonances which decay into D*, are suppressed. 

Before proceeding to a more detailed test, we note that some evidence for the validity 
of the SV limit can be obtained by comparing the inclusive branching ratios measured by 
CLEO 

fir(5 ^ L)*°A) = (0.247 ±0.028) , Br{B ^ D*^X) = {0.239 ± 0.019) . (29) 

These numbers are summed over both charged and neutral B decays, and include semilep- 
tonic and cdu, ccs final states. Both semileptonic and b —>■ ccs decay mechanisms are isospin 
symmetric, so they can not produce an asymmetry between the rates (^). (Although this 
is manifest from (pOD, it holds also for the charge-averaged rates, with multiplicity factors 
added to account for identical particles in the final state.) The only source for such an 
asymmetry is the b — >■ cdu decay mechanism. The large predictions (p!T|) together with 
m imply 

T{Bd D*+X-) + T{B- D*+X—) (30) 
= E |Fo(X)|2 + 3|Fi(X)p + V2Re{Fo{X)F*{X)) + {\h,\' + \g,\') + 0(l/iV,) 

X i=0,l 

r{Bd D*°A°) + T{B- D*^X-) (31) 
= E \F,{X)\^ + 3|Fi(A)p - v^Re(Fo(X)F;(X)) + ^ (j/^.p + \g^\^) + 0{l/K) . 

X i=0,l 

Neglecting the l/N^ corrections, the approximate equality of the measured rates (|9]) can be 
interpreted as evidence for a small ratio Fi/Fq, which coincides with the expectation from 
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the SV limit^. Therefore we will neglect Fi and compare the measured decay spectra against 
the theoretical prediction from factorization for Fq. 

We note a subtlety in applying the combined large Nc and SV limits to inclusive decays. 



It is well known [0 that in the large Nc limit widths are dominated by states with the 
minimum possible number of final state mesons. Each additional meson comes at the price 
of a l/Nc suppression factor. Thus, naively one should include in Eq. ( P^D only states in 
which X is a single meson. However, in the SV limit the energy available for the decay 
becomes arbitrarily large and phase space effects invaliate the naive conclusion. To see this, 
consider large but fixed b and c masses, and let n be the number of single resonances X 
kinematicaly allowed in the decay, that is, lighter than about rrib — rric. The number of 
combinations of m mesons in the final state scales as n"™, where a is a fixed constant^. So 
we see that, in relation to the width into one light meson, the width into m light mesons 
roughly scales as (n/Nc)"^. Thus, there is no suppression of multimeson states for — rric 
large enough so that, roughly, n > N^.- 

Hence in the combined large N^. and SV limits the ratio of the inclusive and semileptonic 
B D*ii' rates can be predicted in terms of the spectral functions Vi and ai for the vector 
and axial currents. The spectral functions are defined through 

{q,qu - gV)n5f(g') + g^g.n^f = %^d'x e''^o\TJl{x).Mo) |o) , (32) 

vi{q^) = 27rImHf=^(g2) , ai{q^) = 27rImHf=^(g2) , 

and have been measured |^ by the ALEPH collaboration in r-decays up to = 3.0 GeV^. 
Above the resonant region the Vi + ai data displays a plateau, in excellent agreement with 
the operator product expansion prediction. Therefore one can safely extrapolate the data for 
> 3.0 GeV^ using perturbative QCD, (t>i + ai) = 1.1. This extrapolation is independent 
of the mechanism responsible for factorization and therefore does not bias our analysis in 
any way. 

For hadronic states X^, not including charmed hadrons the prediction from factorization 

is 



= 3C2|K 



dg^ _ 



— r(5 ^ D*iu) 

dg^ 



^i(g') +01(9') +••• , (33) 



where the ellipses denote terms suppressed by l/N^ or two powers of the SV expansion 
parameter. Note that although the numerical results for B D*ev in the SV limit are 
distorted by the modification in the kinematical factors, such modifications occur in a similar 
way for the inclusive and exclusive decays and therefore cancel out in the ratio. We emphasize 
that Eq. ( p3D can be used to give a very clean factorization prediction for the dr(i? — >■ 
D*Xu)/dq'^ decay rate with input from the r-spectral functions^ and from the measured 
B D*C,v form factor. 



^ We neglected for this argument the non-prompt D* production, which is itself justified in the SV limit. 

Therefore the equality of the rates (29) provides only a consistency check on the validity of the SV limit. 
^ For example, a = 1/2 in the 't Hooft model in the case of decays of a heavy meson into light mesons 19, 20[. 
^ The T decay data in this context was also used in |24| for a test of duality in nonlcptonic B decays. 
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FIG. 1: Inclusive differential decay rate dr(i? — > D^*Xu)/dx. The dashed line shows the result 
in the B rest frame. Here x is the rescaled D* momentum, x = |p/)|/(4.95 GeV) and we normalize 
the rate as in To compare to the CLEO [^] data we boost to the T rest frame which gives 
the solid line. Only three data points (shown) are available that are not contaminated by charm 
contributions above the Dg threshold. 
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FIG. 2: a) Breakdown of contributions to the B — > D*X spectrum in the B rest frame from the 
r (dotted line), sc (dot-dashed line), {e,//} (thin solid line), and du (dashed line) final states. The 
thick solid line shows the total result, b) Inclusive differential decay rate dT/dp for the process 
B D^*X using factorization from the large Nc and SV limits. The two curves are as in Fig. |^. 
Charmed states in X are included using perturbative QCD to "model" the contributions of the 
{sc)v-A spectral density. For large q'^ (small x) this approximation should be reasonable. 



The CLEO data is presented as a function of the rescaled lab momentum of the D* meson 
X = |p|/(4.95 GeV). In Fig. ^ we show the boosted^ spectrum assuming that the B mesons 



^ In Ref. [p5l an analysis of the CLEO data was performed using factorization as in Eq. (0) plus a model 
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are produced mono chromatically and isotropically, with velocity (3 = 0.06 in the lab frame, 
and using 

driab 1 \p\ , 1 dlcM 

ds— -— — . (34) 



d\p\ 2/572 E Js^,„ Pis) ds 



Here 7 = 1/^/1 — (3'^ and E = _|_ mjj,, and P{s) and dFcM/ds are the momentum of the 
D* and the decay rate in the rest frame of the B meson, respectively. The limits of integration 
are Smax = + "^ii. - 2mB'y{E - (3\p\) and Smin = max(0, m| + m|,. - 2mB'y{E + l3\p\)). 
Fig. ^ shows our result for the x spectrum as a continuous curve. The CLEO data is shown 
as data points. The semileptonic rate was taken from the BELLE fit to a form factor 
using a unitarity constrained parametrization P3] . The available data is in good agreement 



with the factorization prediction. 

Unfortunately, for large it is necessary to include charmed hadrons in X to enable a 
comparison with the experimental data that is currently available from CLEO |]TB|. If we 



include charmed hadrons in the X produced through a {sc)v-a current then it becomes 
harder to test factorization in a clean way. One problem is that no data for the {sc)v-a 
spectral functions is available. To make predictions for the {sc)v-a current we model the 
spectral functions using perturbative QCD. This neglects the resonant structures in the 
region above ~ ""^d^) but should be fine for testing factorization at large invariant masses. 
A second potential problem is that the (cs) and (rUr) contributions require b ^ c form factors 
which are not accessible in B ^ D^*^ei' decay. However, since the SV limit implies heavy 
quark symmetry, this can be used to predict these form factors. Finally, an additional source 
of uncertainty is introduced by the value of the charm quark mass. We choose rric = 1.5 GeV 
in our calculations. Note that since the {sc)v-a current is phase space suppressed it only 
contributes roughly half as much to the decay rate relative to {du)v-A current. This helps 
to reduce the model dependence of our B D*X predictions. 

In Fig. we show a breakdown of the lepton and quark contributions to the total B — > 
D*X rate in the B rest frame. The prediction from factorization is in moderate agreement 
with the data. The disagreement does not seem to scale with the invariant mass of the 
system, however due to the large theoretical uncertainties in our calculation, the results in 
the region above the ccs threshold are inconclusive. Improved measurements of the energy 
spectrum in Fig. |I] which disentangle the states with different charm quark numbers would 
help to clarify this issue. 



IV. DISCUSSION AND CONCLUSIONS 

Data suggests that the amplitude for B^ D^n^ factorizes. This can be understood 
via large Nc counting or via large Q as an expansion in powers of 1/Q, the inverse of the 
energy released. Thus, deviations from factorization in this process are doubly suppressed. 
However, data also suggests factorization in the exclusive processes B D*'k^tt~'k~tt'^ and 
B D*uJ7T~ as well as in the inclusive B D*X. Neither large A^^^ nor large Q explain 
by themselves these results. Indeed, large Nc predicts instead a sum of factorizable terms. 
On the other hand the corrections to factorization in large Q are order mx/Q, and in these 
decays, particularly in B ^ D*X, this ratio can be comparable with unity. 



for the X' contributions. Our approach differs from this by the use of the SV limit, r-decay data, and 
the inclusion of the boost to the lab frame. 
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We have pointed out that the condition for large energy release in D^7r~ can be 

understood as a consequence of the SV limit, mb,mc ^ rrib — rric ^ Is.. In the combined 
large Nc and SV limits we can justify factorization m B ^ D*X out to arbitrary invariant 
hadronic mass. 

We suggest the following physical picture. The corrections to factorization are paramet- 
rically small, of order l/N^, but depend on kinematic variables in a way that can amplify 
the magnitude of the corrections if the recoiling particle is not moving fast. The SV limit 
ensures that this kinematic enhancement is absent: for cases for which we have control like 
B^ D^n^, the correction term is ~ l/Q so it depends on the kinematics but in such a 
way as to further suppress the correction to factorization. 

This picture is additionally supported by D and K decay data. When two body D decay 
amplitudes are written in terms of weak transitions into definite isospin states, Aj, and final 
state interaction phases, exp{i6j), it is found that the amplitudes Aj do factorize, provided 
one uses a modified Nc counting Ci(mc) ~ C2{mc) ~ 0(1) [Q. The strong phases, taken 
from experiment, are not small but should vanish in the large limit. We interpret this 
as the expected kinematic enhancement of the correction to factorization, and it suggests 
that factorization fails precisely because the kinematic enhancement shows up mostly in 
6i ~ 1/Nc rather than in Aj ~ {l/Ncf. 

Is this kinematic suppression due to large velocity of the products or large energy of 
the products? In the first case, which corresponds to the color transparency argument, the 
suppression is expected to behave as I/7. In the second the suppression should be given by 
the larger of A/Q or mx/Q, where A is a typical hadronic scale, mx is the invariant mass 
of state produced from the current by the factorizing current, and Q is its energy. We see 
that in K ^ vrvr decays one is bound to have A/Q ~ 1, while 1/7 = 2mj^/mK ~ 1/2 and is 
parametrically suppressed in the chiral limit. Unfortunately, the actual pion mass is too far 
from the chiral limit to distinguish between the two alternatives. 

We studied in this paper the factorization predictions for B D^*^X, focusing on meth- 
ods which can distinguish among the various possible explanations for factorization. In 
Sec. II we derived large Nc relations among isospin amplitudes, which lead to observable 
predictions among decay rates. Imposing additionally the 5*^ limit gives even more rela- 
tions, since certain amplitudes allowed in the pure Nc limit are suppressed. The large Q limit 
by itself gives predictions similar to the combined limit of large A^^^ and SV, however these 
predictions should fail outside of a limited kinematical range. We presented several such pre- 
dictions that can be tested experimentally. Using available data from BELLE and Babar we 
discussed such predictions for the modes B — >■ D^*>Tnr, B D^*^KK and B D^*W-*^K. 

In Sec. Ill we calculated the inclusive decay rate B — > D*X in the combined limit of 
large Nc and SV. In this limit, factorization is expected to work out to arbitrary invariant 
hadronic mass, in contrast to the predictions from large Q factorization in which factorization 
breaking corrections should scale as rrix/Q. The cleanest theoretical prediction is for the 
B — > D*Xu decay, for which we can use r-decay data to extract the spectral function in 
the resonance region^. The comparison of the prediction with the CLEO data is shown in 
Fig. |I| and shows good agreement. Inclusive data for this process is available over a much 



® We take this opportunity to comment on the recent measurement by BELLE ||T^ of the ratio Br{B 
Df'K-)/BR{B- DItt-) = 0.89±0.14, which has been interpreted as a test of factorization. However, 
as pointed out in Rcf. this ratio depends sensitively on unknown subleading Isgur-Wise functions, 
which can accomodate values in the range 0—1.5 within factorization. 
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larger kinematic range, however the data then includes contributions from the {sc)v-a part 
of the current. No data is available for the spectral function of this current, and we therefore 
have to rely on perturbative calculations. The result is shown in Fig. |^, however the large 
theoretical uncertainties preclude us from drawing definite conclusions. 

We emphasize that it would be desirable to separate the b cdu and h ccs final states 
in the B — >■ D^*^X data, including separate measurements for B and D charge eigenstates. 
The methods described in this paper would then allow a more detailed study of factoriza- 
tion in this decay, and should help shed light on the mechanism underlying the observed 
factorization in B decays. 



Acknowledgments 

We are grateful to Vivek Sharma for many useful discussions. This work is supported 
in part by the Department of Energy under contracts No. DOE-FG03-97ER40546 and DE- 
FG03-00-ER-41132. 



[1] C. W. Bauer, D. Pirjol and I. W. Stewart, Phys. Rev. Lett. 87, 201806 (2001) 

[2] C. W. Bauer, S. Fleming and M. Luke, Phys. Rev. D 63, 014006 (2001); C. W. Bauer, 
S. Fleming, D. Pirjol and I. W. Stewart, Phys. Rev. D 63, 114020 (2001); C. W. Bauer and 
L W. Stewart, Phys. Lett. B 516, 134 (2001); C. W. Bauer, D. Pirjol and I. W. Stewart, Phys. 
Rev. D 65, 054022 (2002). 

[3] M. J. Dugan and B. Grinstein, Phys. Lett. B 255, 583 (1991). 

[4] M. Bauer, B. Stech and M. Wirbel, Z. Phys. C 34, 103 (1987). 

[5] M. Neubert and B. Stech, in Heavy Flavors, 2nd ed., edited by A. J. Buras and M. Lindner 

(World Scientific, Singapore, 1998), p. 294. 
[6] H. D. Politzer and M. B. Wise, Phys. Lett. B257, 399 (1991); M. Beneke, G. Buchalla, M. 

Neubert and C. T. Sachrajda, Phys. Rev. Lett. 83, 1914 (1999); Nucl. Phys. B 591, 313 

(2000). 

[7] A. J. Buras and L. Silvestrini, Nucl. Phys. B548, 293 (1999). 
[8] M. A. Shifman, Nucl. Phys. B 388, 346 (1992). 

[9] Z. Ligeti, M. E. Luke and M. B. Wise, Phys. Lett. B 507, 142 (2001). 
[10] K. W. Edwards et al, (CLEO Collaboration), Phys. Rev. D65, 012002 (2002). 
[11] M. A. Shifman and M. B. Voloshin, Sov. J. Nucl. Phys. 47, 511 (1988) [Yad. Fiz. 47, 801 

(1988)]. 

[12] C. G. Boyd, B. Grinstein and A. V. Manohar, Phys. Rev. D 54, 2081 (1996). 

[13] L. Gibbons et al. (CLEO Collaboration), Phys. Rev. D 56, 3783 (1997). 

[14] K. Abe et al, (BELLE Collaboration), BELLE-CONF-0233, BELLE-CONF-0235, presented 

at the 31*** International Conference on High Energy Physics, Amsterdam, Netherlands. 
[15] M. S. Alam et al. [CLEO Collaboration], Phys. Rev. D 50, 43 (1994) [|arXiv:hep-ph/9403295 |. 
[16] BELLE Collaboration, |hep-ex/020704l| . 

[17] BABAR Collaboration, contribution to the 31*** International Conference on High Energy 



Physics, Amsterdam, Netherlands, hep-ex/0207086 



[18] S. Coleman, Aspects of symmetry : selected Erice lectures, Cambridge Univ. Press, 1985 
[19] B. Grinstein and R. F. Lebed, Phys. Rev. D 59, 054022 (1999). 



14 



[20] B. Blok, M. A. Shifman and D. X. Zhang, Phys. Rev. D 57, 2691 (1998) [Erratum-ibid. D 59, 

019901 (1999)]. 

[21] R. Barate et al. (ALEPH Collaboration), Eur. Phys. J. C 4, 409 (1998). 

[22] Z. Ligeti and A. V. Manohar, Phys. Lett. B433, 396 (1998). 

[23] K. Abe et al. (Belle Collaboration), Phys. Lett. B526, 247 (2002). 

[24] C. G. Boyd, B. Grinstcin and R. F. Lebed, Nucl. Phys. B 461, 493 (1996). 

[25] X. Calmct, T. Mannel and I. Schwarze, Phys. Rev. D61, 114004 (2000). 

[26] A. K. Leibovich, Z. Ligeti, I. W. Stewart and M. B. Wise, Phys. Rev. D 57, 308 (1998). 



15 



